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ELECTRON CAPTURE BY PROTONS IN HYDROGEN 
IN AN ELECTRIC FIELD 


SUMMARY 

34324 

The evaluation of the capture cross section in the Brinkman Kramers ap- 
proximation with the hydrogen atom wave function in parabolic coordinates is 
much simpler than this wave function in spherical coordinates. A simple closed 
formula for the electron capture cross section in parabolic coordinates, appro- 
priate to the charge transfer in an electric field, is derived with the electron 
initially and finally in two arbitrary sublevels of the atom, and it is shown that 
the cross section, apart from a weighting factor, is independent of the particular 
sublevel and depends only on the principal quantum numbers involved. In this 
way the charge transfer cross section in an electric field is obtained. In addi- 
tion, the solution of the problem in parabolic coordinates leads, through a trans- 
formation, to the determination in closed form of the cross section is spherical 
coordinates. A general formula is worked out for the cross section o-(n/,n'/')> 
and as an example simple formulas for cr (Is, in) with the final states the sub- 
levels of n = 1, 2, 3 are tabulated. Previous works have listed these cross sec- 
tions in integral forms. 

It is found that the momentum distribution of electrons in the Stark levels 
of the hydrogen have extremely simple form in contrast to the momentum dis- 
tribution in the optical levels. The averaged momentum distribution over the 
sublevels of a principal quantum number gives results identical to the results 
originally obtained by Fock, and recently reobtained by a different method by May. 

The classical momentum distribution of a bound particle in a Coulomb field 
is shown to be the same in form as the quantum mechanical distribution. The 
use of this distribution function in classical calculation of scattering will make 
the classical and the quantum mechanical calculations more consistent. * il 
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I. INTRODUCTION 


The momentum density distribution function for a bound electron in the 
Coulomb field of a nucleus of charge Ze occupying the state of the principal 
quantum number n and the azimuthal quantum number / has been found by 
Podolanski and Pauling 1 using the Fourier transform of the spatial wave functions, 
and by Fock 2 by solving the Schroedinger integral equation in the momentum space. 
Fock further has shown that when the above density distribution function is 
averaged over l , a simple function results. Let p represent the momentum of 
the electron and q its propagation vector, then p ="hq and this function is given 
by 3 


D(n, q) = 


n 


n - 1 

L 

1=0 ■ 


(2/ + 1) D(n/ , q) = 


8 Z 2 a 3 

n 2 77 2 a 2 


1 


2x4 


O + q ) 


a = — ( 1 ) 
na 0 

with a Q the Bohr radius. Eq. (1) has recently been rederived by May 4 by a 
different method. 

The momentum density distribution function of a bound electron in a Coulomb 
field has a particular simple form when the wave function is expressed in para- 
bolic coordinates, which is the wave function of the hydrogen atom in an electric 
field, weak compared to the electron-nucleus interaction. If we specify a particular 
state of the electron bynn n 2 , with n x and n 2 the quantum numbers corresponding 
to the i and 17 coordinates, it is shown here that the momentum distribution 
function, except for a weighting factor, has the same form as D(n, q) given by 
Eq. ( 1 ), while the same function describing a state n/ in spherical coordinates 
has a more complicated form. It is given in terms of the Gegenbauer polynomials 
whose argument is a function of the q 5 . 
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The electron capture cross section by proton from the hydrogen atom in 
the Born approximation in which the interaction between the nuclei is neglected 
is given as an integral over the product of the squared moduli of the Fourier 
transforms of the initial and the final electron states wave functions. 6,7 This 
integral is the same as an integral over the product of the momentum densities 
of the two states. This is the approximation used in this paper and by using 
momentum density functions of the Stark levels of the hydrogen atom which are 
similar to the density function given in Eq. (1), the integral is found in closed 
form. Furthermore it is found the interesting result that the capture cross 
section for the initial state nrijn 2 and the final state n'nj n' is independent of 
n r r» 2 , n'j and n' , and depends only on the principal quantum numbers n andn'. 
The other vairable parameters of the cross sections are the charges of the 
target nucleus and the incident particle, Z and Z' , and the initial relative velocity 
of the two heavy particles. 

In problems in which radiative transitions affect the population densities of 
the excited states of the hydrogen atom, and also in the capture problems in the 
non-hydrogenic atoms with wave functions which can be approximated by the 
products of the one electron wave functions, it is desirable to have the capture 
cross sections in spherical coordinates. Bates and Dalgarno 8 have found these 
cross sections between the ground state as the initial state and the sublevels of 
the n = 1 - 4 as the final states. Hiskes 9 has extended these calculations to sub- 
levels through n = 15 as the final states. In all these calculations, however, the 
results are in integral form. By transforming the hydrogenic atom wave func- 
tions from the parabolic to the spherical coordinates wave functions, and using 
the electron capture amplitudes in parabooic coordinates, it is possible to express 
the cross section in spherical coordinates in closed form, and as a sum over 
finite number of terms. This has been done here and a general expression for 
the cross section for the initial state n l and the final state n ' l ' is derived. 

In particular the cross section between the ground state and all sublevels 
through n = 3 are listed. 

At the end, the classical momentum distribution of a bound particle in a 
Coulomb field is worked out in Appendix I, and it is shown that the distribution 
function is similar to the quantum mechanical distribution function given by (1). 

By applying the Bohr-Sommerfeld quantization rule to the classical function, 
the two functions become identical. The consistency of the classical and the 
quantum mechanical distribution functions makes it desirable that in a classical 
calculation of some scattering problems the velocity distribution function given 
in Appendix I be used, and comparison be made with the available quantum 
mechanical calculations. 
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II. MOMENTUM DISTRIBUTION IN THE STARK LEVELS 
OF THE HYDROGEN ATOM 

The Schroedinger equation for a bound electron with position vector r in the 
field of a charge Ze located at the origin is given by 


1 z 1 i 

V 2 - i//(a, r) = a 2 i/< a, r), 

2 a o r -l 2 


( 2 ) 


with a defined in (1). The Fourier transform of 0(a, r) is given by 


u(a, q) = (27r)' 3/2 


j* e l9 ' r i//(a,r)dr. 


(3) 


We define in addition the function 


v (a, 


q) = (2 t r) 


- 3/2 


Je ia ' r | '/'(a, 


r) dr. 


(4) 


Then by substituting (Z/r)<//(a, r) from Eq. (2) into Eq. (4) and carrying out a 
partial integration we obtain 


«(a, q) 


2 

a 0 ( a2 + d 2 ) 


y (a, q). 


(5) 


The momentum density, and also the electron capture cross sections, are 
given in terms of the squared modulus of « (a, q) summed over the magnetic 
quantum numbers. Since this quantity is a scalar and invariant under rotation 
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of the coordinate system, for evaluation of «( a, q) we can take for convenience 
the z-axis of the coordinate system along the q and we designate in this case 
the «(a, q) and v(a, q) by u(a , q) and v(a, q) respectively. 

We now evaluate v(a, q) in the parabolic coordinates. Keeping in mind that 
in these coordinates r = 1/2 + 17), z = 1 / 2 - 17), dr = 1/4 (£ + 77) dg dr] d 4 >, 

and expressing the spacial wave function as the product of two associated Laguerre 
functions 10 , we obtain from ( 4 ): 


i^nn^m, q) = (27 r)‘ 3/2 j*e iqz — P (nnjn 2 m, r) dr = 0, m/ 0 , 


( 6 ) 


z^nrijnjO, q) = ( 277 ) 


-3/2 


f-i 


i/'(nn 1 n 2 o, r) d] 


1 

Za 3/2 

/* 00 /*00 

exp 

•'o •'O 

r a i Q 

--(^ + 7?)+— (£-*?) 

2 2 J 

L n,( a O L° 2 (a7 ] )d^d77 

277 /2n 

n i !n 2 ! 

1 

Za 3/2 

I( ni ) I*(n 2 ), 


(?) 


277/2n rij ! n 2 ! 


where 


1(0 



exp 




L° (aO d£. 


( 8 ) 


Using the generating function of the associated Laguerre functions 11 it follows 
from (8) that 


4 



s' 1(0 _ 1 


dx 


i 


00 




-1 


— (a - iq) (1 - s) + as 
2 


( 9 ) 


Now by introducing 


co = — (a - iq), 
2 


( 10 ) 


the preceding equation can be written 


E 


8 * 1(0 = 1 
/ | O) + <^>*s 



( 11 ) 


By equating the coefficients of equal powers of s on both sides of this equation 
we obtain 


1(0 


jM 

CO 



( 12 ) 
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Substitution of !(/) from this equation into Eq. (7) gives 


t^nn^o, q) = 

27r/^n 


Za 3/2 


co 



( 13 ) 


By making use of the equation 


n i + n 2 = n - 1 - m, 


(14) 


Eq. (13) can be written, after neglecting a constant phase factor, 


v(nn 1 n 2 o, q) 




(15) 


and by Eqs. (5, 10), 


a(nn 1 n 2 o, q) 


1 Za 3/2 

277 /2TT 2 a 0 M 4 



(16) 


Equations (15, 16) are the main equations from which the momentum distribution 
and electron capture cross sections are derived. 
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The momentum distribution of an electron in the state nri^ is given by 


DtnnjOj, q) ■ = | s^nn^m, 


q) 


“(nrij^o, q)| 2 


8 Z 2 a 3 
n 7 t 2 a 2 


1 

(a 2 + q 2 ) 4 ’ 


(17) 


an expression which is independent of the sublevel quantum numbers n x and n 2 . 


The averaged momentum distribution of an electron in the shell n is given 
by averaging D(nn x n 2 , q) over all n x and n 2 . Keeping in mind that there are n 2 
states for a given n we obtain 


D(n, q) = n 


n- 1 


L 

n i =n 


D(nnjn 2 , q) = n' 1 D(nn x n 2 , q) 


8 Z 2 a 3 1 

= • . ( 

n 2 tt 2 a 2 (a 2 + q 2 ) 4 

It is interesting to notice that the integral with respect to q over D(n, q) is 
normalized to unity which corresponds to an electron in the n shell: 


f D(n, q) dq = f 

J n 2 v a 2 J 0 


q 2 dq 
(a 2 + q 2 ) 4 


32 Z 2 a 3 . 77 _ 

n 2 77 a^ 32 a 5 


(19) 


The expression for D(n, q) given in Eq. (18) is identical to the expression 
given in Eq. (1), derived by Fock and May using spherical coordinates, and to the 
expression given in Eq. (All) for the classical momentum distribution in a 
Coulomb field. 
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Equation (17) has the interesting meaning that the Stark levels of the hydrogen 
atom within a given shell n, except for a weighting factor, have the same momen- 
tum distribution. 


III. CAPTURE CROSS SECTION IN AN ELECTRIC FIELD 12 

Assume a nucleus of charge Z ' and mass M' captures an electron from a 
single electron atom with nucleus of charge Z and mass M to form an atom with 
nucleus of charge Z' and mass M'. The capture cross section in the Born ap- 
proximation in which the interaction between the nuclei is neglected is given 
by 8 - 1 3 


'! i. /) 


2tt /j. 2 e 2 

h 4 


-f 

v J-i 


|f(a, q) | 2 x |g(a\ q ')| 2 d(cos 6), ( 20 ) 


with p. the reduced mass of the system, e the absolute value of the electronic 
charge, v and v' the magnitudes of the initial and final velocities of relative 
motion of the nuclei; and 6 the angle between v and v ' . From here on the un- 
primed symbols correspond to the initial states while the primed are for the 
final states. The functions f( a , q) and g( a , q) are defined by 


f (a, q) 


g(a, q) 


J 

J 


e iQ-r 


e id-r 


i//(a, r) dr, 


— '/'(a, r) dr, 

r 


( 21 ) 


( 22 ) 


with \p(a, r) and i p(a , r") the wave functions of the atoms Z and Z' , and q 
and q' are related to the velocities v and v' through 


with 


q = 


M 

M + m 


k - k' 


q 


= k - 


M' 

M' + m 


k', 


M'(M + m) 
ti(M' + M + ml V ’ 


M(M' + m) 
ti(M' + M + m) 


f 

v . 


(23) 


(24) 
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In the impact parameter approximation the corresponding expression for 
the cross section is given by 4, 14,15 

V 


cr(i, f) 



f(a, q) | 2 x !g(a', q')| 2 pdp. 


with p a dimensionless quantity defined by 


(25) 


p = tiv/ e 2 , 


and 


with 


q 2 = p 2 + q 2 > q ' 2 = p 2 + q ' 


• 2 


q, = P = — 

2a 0 P 


p - 


' z 2 Z ' 2 


, n" n 


#2 


q! 


-1 


2a oP 


fz 2 Z ' 2 

p 2 + b - ^ 


(26) 


(27) 


(28) 


with a Q the Bohr radius and n and n' the principal quantum numbers of the 
initial and the final states of the electron. 


As the integral that appears in the impact parameter method is easier to 
handle, the cross sections evaluated below are according to this method. 

In a weak electric field the energy levels of the hydrogen atom split into 
Stark levels. The corresponding wave functions are given in parabolic coordinates. 
By the definitions (21, 22) and through Eqs. (15, 16) we obtain 


f(nn x n 2 o, q) 


g(nnjn 2 o, q) 



Za 3/2 /q ) *\2n 1 

2 a 0 |co| 4 U/ 


Za 3 / 2 


&T' 


(29) 


(30) 


The capture cross section in the impact parameter approximation with the atom 
Z in the state nri^ and the atom z ' in the state n'nj n 2 is then given by 
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crCnrijn^n'njn') 


27tZ 2 Z <2 (ga*) 3 f°° p dp 

16 nn' p 2 Jq M 8 x 1^*1 4 


Since through Eqs. (10, 1, 27, 28) we have 


\co 


'I 2 = i (a' 2 + q' 2 ) = \ (a 2 + q 2 ) = \co\ : 

4 


(31) 




(32) 


it follows that 


crfnnjrij, n'rij'n') = 


2 9 tt (ZZ 1 ) 2 (aa') 3 


nn 


f 

•'o 


pdp 


[a 2 + /3 2 + /? 2 ] 6 


— 2 2‘ (Z 2* ) s 1 

5(nn') 4 p 2 [(a 2 + /3 2 ) a 2 ] 5 


(33) 


It is evident that as stated previously cr(nn, n 2 , n'nj n') is independent of njn 2 njn* . 
The transition o-(nn 1 n 2 , n’) is obtained by summing cfnn^j, n'njn') overran'. 
Noticing that only states with m* = 0 contribute to the cross section we obtain 


cr(nn 1 n 2 , n') = n' ^nn^j, n'njn') 


= „.» ■ 1 ( 34 ) 

5 n 4 n' 3 p 2 [(a 2 +q 2 )a 2 ] 5 

Finally the transition o-(n, n') is obtained by averaging cr(nn,n 2 , n') over the 
initial states. Since there are n 2 initial states but only n states have nonvanish- 
ing cross sections we obtain 


cr(n, n') = n _1 afnrijnj, n') 


t r a 2 2 8 (ZZ*) 5 . 


5 n 5 n' 3 p 2 [(a 2 +q 2 )a 2 ] s 


( 35 ) 


The above result is identical with the result obtained by May. 
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IV. Capture Cross Section in Spherical Coordinates 

Let cp(n / m | r) represent the wave function of the hydrogen atom in spherical 
coordinates. Then we define the functions F(n l m, q) and G(n l m, q)by 


F(n l m, q) 


G(n / m, q) 



9(n l m, r) dr , 


cp(n l m, r) dr. 


(36) 


(37) 


By writing 


cp(n/ m, r) = R n/ (r) (?), (38) 

and using the expansion 16 

00 k 

exp [i q*r] = 4 n LL ik ’k (qO (i) y^ k) (39) 

k = 0 k 

where r and q are unit vectors in the directions of q and r, Eq. (36) can be 
written as 

F(n/ m, q) = 47Ti / I n/ (q) Y /m (q) , (40) 

where 

i n /(q) = f R„/ j i < qr > r2 dr (41) 

Jo 

When q is along the z axis, by the definition of Y / (q) in terms of the associated 
Legendre polynomials we have 


Y <q = £) = JlLll S(m, o). (42) 

/m ! 4n 
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If in this case as usual we designate the F(n l m , q) by F(n l m , q), through 
Eq. (40) we obtain 


i PTi + 1 

F(n Z m , q) = S (m, o) 4 tt i y — I n ; (q), 


and 


(43) 


F(n l m, q) 


/ 4w 

V 21 + 1 


F(n / o, q) Y u (q). 


When Eq. (44) is summed over m and use is made of 

T |y / m (i)|2 = — 

L I 1 m 477 

m =-/ 

we obtain 


17 


/ 

L 

m=- l 


F(n l m, q) | 2 = F(n l o, q). 


(44) 


(45) 


(46) 


This was the result stated at the beginning of Sec. II by the statement that 
the left hand side of Eq. (46) is a scalar and invariant under rotation of the co- 
ordinate system. 

Since the wave function of the hydrogen atom forms a complete set in each 
of the two spherical and parabolic coordinates, the wave function in spherical 
coordinates can be expressed as a linear superposition of the wave functions in 
parabolic coordinates. In this transformation the principal and the magnetic 
quantum numbers n and m remain fixed: 


cp(n l m, r) = 


n - 1 - m 



a / ni '/'(nn,n 2 m ) r), 


(47) 


where a /n are the elements of the transformation matrix. In our problem m = 0. 
By dropping nm from the indices in Eq. (47), replacing by i , substituting 
this equation in Eq. (36) and making use of Eq. (21) we obtain 


n-l 

F (/, q) = 

i = 0 



f(i, q)- 


(48) 
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The determination of for m = 0 and arbitrary n is worked out in Ap- 
pendix II, and the matrices for n = 1, 2, 3, 4, 5 are given explicitly. 

By Eq. (29) we now obtain 


n- 1 n-1 


If (z, q) 1 2 = ZZ a/i a /j f (i> q) f *(i> q) 


i=0 j =0 


77 Z 2 a 3 
4a 2 0 n M 8 


Z */i a <, 


i * J 


co 


* 4 I i - i 


M 


4 i-j 


7T Z 2 a 3 


4na o i^i 8 i,j 


Z a,: a,i 


R 




CO 


4| i- j 


(49) 


where R stands for the real part of 


* 4 | i - j 


When a binomial expansion is made of this quantity in terms of a and q , 
(Cf. Eq. (10)), we obtain 


n-l n-l 2 I i - j 


iki.di»= 64 7,2 ° rr 

naJ(“ 2 + <I J ) 4 4^*4 ^ 

In similar way, through Eq. (30), 




2k 


(a 2 + q 2 >2 | i - j | 
(50) 


|G(/, q) 1 2 = 


_ 16 77 Z 2 a 3 V 1 


,2 , \2 L l 


(~) k a u a n 


^4 1 i - i | \ a 4|i - j| ' 2X q 2 ^ 


2 ^ n2 \2 I i - j 


(51) 


n(a*+q z ) z ' \ 2W (a 2 +q 2 ) 

By making use of Eqs. (27, 28) and making the expansion 


q 2 ^ = ( p 2 + ,6 2 )^ = ^ (^) /S2( ^ ) ^ (52) 

/x.=0 ' ' 

in Eqs. (50, 51), keeping in mind Eq. (32), we obtain the capture cross section 
with the initial states n l and the final states n' / ' , averaged over the initial 
states and summed over the final states. This is given by 
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o-(n/,n'/')= 1 |F(/, q)| 2 x |G(/', q' ) | 2 pdp 

277(2/ + l)p 2 J 0 


pi (M + M ) + l dp' 


2 9 77(ZZ*) 2 (a a') 3 V V" B(y)B(r') f 

(2/ + 1) a 2 p 2 nn' y y J o . p2 ao 2^(3 + I - J I + I i'- i ' I ) 


(a 2 +/3 2 +p 2 ) 


where 


(53) 


D , u f sK /4 1 i - j I 4|i-j|-2\ 2 (X-M) 

w = (-) II J a e 


(54) 


y stands for the set of 4 integers i , j , k, p. which take on the following values: 


i = 0, 1, 2 n - 1 ; 

j = 0, 1, 2, ..... n - 1 ; 

X = 0, 1, 2 2 | i - j | ; 

P = 0, 1, 2 k. 


~s 


> 


J 


The integral appearing in Eq. (53) is of the form 

r°° p 2M+ 1 dp 
I (M, N) = , 

Jo (b 2 + p 2 ) N 

which after integration by parts yields 

/ N- IV 1 

I(M, N) = ( M ) I(0 ’ N “ M) 


(55) 


(56) 


= — (M+l Y 1 f N-1 V 1 b- 2 ( N * M - 1 ). (57) 

2 \ M + 1/ 


14 


In this way we obtain 


cr( n4, n' 4') = tt a* 


2 2 8 (ZZ') 5 


(24 + 1) (nn') 4 p 2 [(a 2 + /3 2 ) a 2 ] 5 


B(y) B(y ' ) 


y y (1 + ix + /x 


) + 2 ( |i - j | + |i' - j ' | )Yy + /S 2 ) 2( | i 

\ 1 + /U. + >Lt / 


i-il + 


By introducing the dimensionless quantity 


C(y) = (a 2 + /3 2 )" 2 l i -il B(y) 


( .,r„ .. 

U \ 2\ / \fi) (a 2 + 


we finally get 


.0 ,0,. 2 2 8 (ZZ*) 5 

cr^n-t., n 4 ) = 77 a„ v 


(24 + 1) (nn') 4 p 2 [(a 2 +/3 2 ) a 2 ] 5 


(1 +M + M) -1 ( 5+2C|l - j| + ^ " J !) j C(y)C(y') 
1 r 1 \ 1 + /j, + fj. / 


y y 


with y defined by Eq. (55), and a and /? given by Eqs. (1) and (16). 

Since cr(n4, n'4') does not change by the interchanges of i and j , and i' 
and j Eq. (60) can be written alternatively as 


/ s 2 2 8 (ZZ') S 

_cr(n4, n 4 ) = na„ - — 


(24 + 1) (nn*) 4 p 2 [(a 2 + /?) a 2 ] ; 


j)] [2 - S(i', j ')] (1+^ + M'r 1 ( 5 + 2C 1 1 + ^ + + \ 3 ^C(y) C(y') , 


y y 
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where S(i, j) is the Kronecker delta, and summation over the indices is the 
same as given by Eq. (55) except that now 


j = 0, 1, 2, • • • • , i; 


j' = 0, 1, 2, 


■ • , i 


(55’) 


If N be the number of the terms in the summations expressed by y in Eq. 
(60'), it is easy to show that 


N = A. n 2 (n + 1) (n + 2). 

D 


(61) 


Then the total number of terms in Eq. (60') is NN\ 

When the initial state is the ground state, Eq. (60') reduces to, omitting the 
primes for convenience, 


cr( 1 s, rrl) - 


it a n 


2 8 (ZZ') S . 1 

n 4 p 2 [(a 2 +/3 2 ) a 2 ] 5 


n" 1 


2(i-i> \ 


LULL « 

i “0 j=0 \= 0 /x=0 


(-> 


)] 


k f4(i - i)' 

2k 


k' 

M / or 


4(i-j )-2k /g 2(k-M) (62) 

(1 + 4 5+2Ci_j) ) (^ + ^) 2(i - i)+ ^ 


1 + /~T 


The following simple formulas are obtained for the electron capture cross 
sections with the ground state as the initial state and the sublevels of n =1, 2, 3 
as the final states. Introducing 


a 2 + /3 2 ’ 

with a and jB referring to the final states (Cf. Eqs. (1, 28), and 

a 0 = 2 8 (Z/Z') s p-2 A 8 (tt a 2 ), 


(63) 


(64) 
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we obtain 


(i) a(ls, Is) =|cr 0 , (65) 

in agreement with the value given by Bates and Dalgamo 8 . 

(ii) o-(ls, 2s) =2 a Q ^|--yA+yA 2 ^> ( 66 ) 

cr( Is, 2p) = 2 a 0 A -yA 2 j, (67) 

cr(l, 2) = cr(l S, 2s) + cr(ls, 2p) = ^ Cr Q . ( 68 ) 


(iii) .(1 s, 3s) = 3a # (| - L6 a + A 2 - ^ A 3 + 256 ^ (69) 

-ds, 3p) = 3a 0 ^ A - A 2 + 32 A 3 - 1|» A 4 ^ , (70) 

a(ls, 3d) =3a 0 -f\* +^A 4 ), (71) 

-(1, 3) = a(l s, 3s) + c( 1 s, 3p) + a( 1 s, 3d) a Q . (72) 
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APPENDIX I 


Classical Momentum Distribution of a Bound Particle 
in a Coulomb Field 

Assume a bound particle of charge -e, mass m and momentum p in the field 
of a nucleus of charge Ze. The conservation of energy requirement is that 


>2 

2m 


e = £1-5s! 


(Al) 


with E the energy of the particle and r its distance from the origin. We recog- 
nize that the density of the particle in phase space is given by 


NS 


(•■fy-r} 


with N the normalization constant determined by the equation 


n / s ( e -§^ +5 ? ! ) drd,, = 1 ' 


(A2) 


On the other hand if the momentum density between p and p + dp is given by 
D(p), we must have 


f- 


D(p) dp = 1 
Combining Eqs. (A2) and (A3) we obtain 


D(p 1 = N| 8 ( E -|i + ^) dr 


(A3) 


j 5 ( E -rl + ^) drd > 


(A4) 
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Below we evaluate the integrals occurring in (A4). Introducing 


„ p 2 Ze 2 

x = E - + 

2m r 


(A5) 


we obtain 


K*-k + ¥) dr = 4 ’ f 


E-p 2 / 2m 


(Ze 2 ) 3 8Q0 
T r. P 2 1 4 

[ x - E + d 


dx 


477 (Ze 2 ) 3 

'eI.eI 4 

_ 2m 


Consequently 


J§ (e -EL + dr dp = 16 77 2 (Ze 2 ) 3 J 


P 2 dp 

,o l ^-*] 4 

L2m J 


(A6) 


= 16 77 2 (Ze 2 ) 3 • 7r(2m) . 3/2 . 

32 (-E) 5/2 


In this way we obtain 


D(p) = 8 (- 2l " E > S/i , 

7r 2 (p 2 - 2m E) 4 


(A7) 


(A8) 


which is the desired relationship. It is seen that D(p) is a scalar and independ- 
ent of the direction of p. 

To show the agreement of this equation with the quantum mechanical momen- 
tum distribution, following the text we introduce the propagation vector q related 
to p by 


p = ti q. 


(A9) 
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We then obtain 


D(q) dq = D(p) dp = ti 3 D(p) dq. 


(A10) 


According to Bohr-Sommerfeld classical quantization formula: 


E 


Z 2 e 2 


n 2 a„ 


ft 2 a 2 
2m ’ 


a defined by Eq. (1). In this way we obtain from (A8), (A9), (A10), 


D(q) 


8 ft 3 (fia) 5 _ 8 Z 2 a 3 

7T 2 "h 8 (a 2 + q 2 ) 4 n 2 7 T 2 a 2 (a 2 + q 2 ) 4 


in agreement with Eq. (1). 


(All) 
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APPENDIX II 


The Elements of the Transformation Matrix 
between the Spherical and the Parabolic Wave Functions 


We require the transformation matrix which for a given principal quantum 
number n transforms the n states with m = 0 in the spherical coordinates to 
another set of n states in the parabolic coordinates. The simplest way to find 
the transformation matrix is by noticing that the matrix elements of interaction 
of an electric field with the hydrogen atom, taken with the unperturbed wave 
function of the atom, are diagonal in the parabolic coordinates representation. 
Let € be this matrix and let A be the transformation matrix. We must then 
have 


AH*A" 1 =€, (A 12) 

where H' is the interaction matrix in spherical coordinates. € is given by 18 


€ ij = “ j) |ne F a 0 (2j + 1 - n), 


j = 0, 1, 2, 


' • • . n 


(A 13) 


with F the strength of the electric field. Similarly through Eq. (38) we can write 


H! . = - eF 


j^iO < r) r COS ^njO^ d] 


CO 

= - eF J R ni< r > R nj (r) r 3 dr x J Y. 0 (Q) Y j0 (fi) cos 0 dfi. (A14) 
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Now 19 


j' v i 0<' n) Y j0< n > cos Sdn-2n (2i±i.) I/J J“ P^P/Wd,. 


1 

2 



[i Pf., p® + (i + 1) p® +1 p«] d M 


iS(i, i + 1 ) + (i + 1 ) 8 (i + 1 , j) 

[( 2 i + 1 ) ( 2 i - 1)] 1/2 [( 2 i + 1 ) ( 2 i + 3 )] " 2 ‘ 


Also 20 



R(. r > Rt.O, 

ni ni"l 


dr=|a 0 n 


/n 2 " 


- l" 


In this way we obtain 


(A 15) 


(A 16) 


H'. = - - ne F a Q [C. S(i, j + 1 ) +C . +1 S(i + 1 , j)], (A 17 ) 


with 


C. 


i / 7 TI 2 


[C2i + 1> C2i - 1>] 1/2 
Eq. (A 12) can now be written 


2>a H ;,~uV 


(A 18) 


(A 19 ) 


which by means of Eqs. (A13), (A17) reduces to 


j+i S + i + j-i Cj - (2i + 1 - n) A . . . 


(A20) 
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This is a recursion relationship for A by means of which all the elements of a 
row of A can be found once the first element of the row is given. 

The unitary condition on A gives 

n~l 

(A21) 

j-0 

where we have assumed that the elements of A are real and if a solution for A 
is found our assumption is justified. 



Eqs. (A20) and (A21) are sufficient for the determination of A for any given 
n. This has been done for n = 1, 2, 3, 4, 5. In the text the inverse of A is needed. 
Let a be the inverse of A, 

a = A -1 = A , (A 22) 

with A the transpose of A . The elements of a for n = 1, 2, 3, 4, 5 are given 
below. 
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